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1 Introduction 



Quantum gravity remains an outstanding problem of fundamental physics. The bottom 
line is we don't even know the nature of the system that should be quantized. The 
spacetime metric may well be just a collective description of some more basic stuff. The 
factP] that the semi-classical Einstein equation can be derived by demanding that the 
first law of thermodynamics hold for local causal horizons, assuming the proportion- 
ality of entropy and area, leads one to suspect that the metric is only meaningful in 
the thermodynamic limit of something else. This led me at first to suggest that the 
metric shouldn't be quantized at all. However I think this is wrong. Condensed matter 
physics abounds with examples of collective modes that become meaningless at short 
length scales, and which are nevertheless accurately treated as quantum fields within 
the appropriate domain. (Consider for example the sound field in a Bose-Einstein con- 
densate of atoms, which loses meaning at scales below the so-called "healing length", 
which is still several orders of magnitude longer than the atomic size of the fundamental 
constituents.) Similarly, there exists a perfectly good perturbative approach to quan- 
tum gravity in the framework of low energy effective field theory i2^. However, this is 
not regarded as a solution to the problem of quantum gravity, since the most pressing 
questions are non-perturbative in nature: the nature and fate of spacetime singularities, 
the fate of Cauchy horizons, the nature of the microstates counted by black hole entropy, 
and the possible unification of gravity with other interactions. 

At a shallower level, the perturbative approach of effective field theory is nevertheless 
relevant both for its indications about the deeper questions and for its application to 
physics phenomena in their own right. It leads in particular to the subject of quantum 
field theory in curved spacetime backgrounds, and the "back-reaction" of the quantum 
fields on such backgrounds. Some of the most prominent of these applications are the 
Hawking radiation by black holes, primordial density perturbations, and early universe 
phase transitions. It also fits into the larger category of quantum field theory (qft) 
in inhomogeneous and/or time-dependent backgrounds of other fields or matter media, 
and is also intimately tied to non-inertial effects in flat space qft such as the Unruh 
effect. The pertubative approach is also known as "semi-classical quantum gravity", 
which refers to the setting where there is a well-defined classical background geometry 
about which the quantum fluctuations are occuring. 

The present notes are an introduction to some of the essentials and phenomena of 
quantum field theory in curved spacetime. Familiarity with quantum mechanics and 
general relativity are assumed. Where computational steps are omitted I expect that 



4 



the reader can fill these in as an exercise. 

Given the importance of the subject, it is curious that there are not very many 
books dedicated to it. The standard reference by Birrell and Davies0 was published 
twenty years ago, and another monograph by Grib, Mamaev, and Mostapanenko|3], half 
of which addresses strong background field effects in fiat spacetime, was published two 
years earlier originally in Russian and then in English ten years ago. Two books with a 
somewhat more limited scope focusing on fundamentals with a mathematically rigorous 
point of view are those by Pulling [S] and Wald[6^. This year DeWitt[7] published a 
comprehensive two volume treatise with a much wider scope but including much material 
on quantum fields in curved spacetime. A number of review articles (see e.g. |H1 El HHl 
nil 1121 IH] ) and many shorter introductory lecture notes (see e.g. [HI [T31 HH] ) are also 
available. For more information on topics not explicitly referenced in the text of these 
notes the above references should be consulted. 

In these notes the units are chosen with c = 1 but h and G are kept explicit. The 

spacetime signature is (H ). Please send corrections if you find things here that are 

wrong. 

2 Planck length and black hole thermodynamics 

Thanks to a scale separation it is useful to distinguish quantum field theory in a curved 
background spacetime (qftcs) from true quantum gravity (qg). Before launching into 
the qftcs formalism, it seems worthwhile to have a quick look at some of the interesting 
issues that qftcs is concerned with. 

2.1 Planck length 

It is usually presumed that the length scale of quantum gravity is the Planck length 
Lp = {hG/c^Y^"^ ^ 10^'^^ cm. The corresponding energy scale is 10^^ GeV. Recent 
"braneworld scenarios" , in which our 4d world is a hypersurface in a higher dimensional 
spacetime, put the scale of quantum gravity much lower, at around a TeV, corresponding 
to I/Tev = lO^^Lp ^ 10~^^ cm. In either case, there is plenty of room for applicability 
of qftcs. (On the other hand, we are much closer to seeing true qg effects in TeV scale 
qg. Por example we might see black hole creation and evaporation in cosmic rays or 
accelerators.) 

Here I will assume Planck scale qg, and look at some dimensional analysis to give a 
feel for the phenomena. Pirst, how should we think of the Planck scale? The Hilbert- 
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Einstein action is She = (^/167rLp) J d'^x\g\^^'^R. For a spacetime region with radius 
of curvature L and 4-volume the action is ~ hlL/Lp)"^. This suggests that quantum 
curvature fluctuations with radius less than the Planck length L < Lp are unsuppressed. 

Another way to view the significance of the Planck length is as the minimum lo- 
calization length Ax, in the sense that if Ax < Lp a black hole swallows the Ax. To 
see this, note that the uncertainty relation AxAp > h/2 implies Ap > h/Ax which 
implies AE > he/ Ax. Associated with this uncertain energy is a Schwarzschild ra- 
dius Rs{Ax) = 2GAM/c^ = 2GAE/c^, hence quantum mechanics and gravity imply 
Rs{Ax) > Lp/Ax. The uncertain Rs[Ax) is less than Ax only if Ax > Lp. 

2.2 Hawking effect 

Before Hawking, spherically symmetric, static black holes were assumed to be completely 
inert. In fact, it seems more natural that they can decay, since there is no conservation 
law preventing that. The decay is quantum mechanical, and thermal: Hawking found 
that a black hole radiates at a temperature proportional to h, Th = (^/27r)K, where k, 
is the surface gravity. The fact that the radiation is thermal is even natural, for what 
else could it hel The very nature of the horizon is a causal barrier to information, and 
yet the Hawking radiation emerges from just outside the horizon. Hence there can be 
no information in the Hawking radiation, save for the mass of the black hole which is 
visible on the outside, so it must be a maximum entropy state, i.e. a thermal state with 
a temperature determined by the black hole mass. 

For a Schwarzschild black hole k = 1/4GM = l/2Rs, so the Hawking temperature 
is inversely proportional to the mass. This implies a thermal wavelength = Sn'^Rg, a 
purely geometrical relationship which indicates two things. First, although the emission 
process involves quantum mechanics, its "kinematics" is somehow classical. Second, 
as long as the Schwarzschild radius is much longer than the Planck length, it should 
be possible to understand the Hawking effect using only qftcs, i.e. semi-classical qg. 
(Actually one must also require that the back-reaction is small, in the sense that the 
change in the Hawking temperature due to the emission of a single Hawking quantum 
is small. This fails to hold for for a very nearly extremal black holejTTj.) 

A Planck mass (~ 10^^ gm) black hole — if it could be treated semi-classically — 
would have a Schwarzschild radius of order the Planck length (~ 10~'^^ cm) and a 
Hawking temperature of order the Planck energy (~ 10^^ GeV). From this the Hawking 
temperatures for other black holes can be found by scaling. A solar mass black hole 
has a Schwarzschild radius i?^ ~ 3 km hence a Hawking temperature ~ 10~^^ times 
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smaller than the Planck energy, i.e. 10 GeV. Evaluated more carefully it works out 
to Th ~ 10^'' K. For a mini black hole of mass M = 10^^ gm one has Rg ~ 10"^^ cm 
and Th ~ 10"K ~ 10 MeV. 

The "back reaction", i.e. the response of the spacetime metric to the Hawking 
process, should be well approximated by the semi-classical Einstein equation Gfj,i, = 
87iG{T^^) provided it is a small effect. To assess the size, we can compare the stress 
tensor to the background curvature near the horizon (not to G^iy, since that vanishes in 
the background). The background Riemann tensor components are ~ 1/i?^ (in suitable 
freely falling reference frame), while for example the energy density is ~ Tfj/h^ ~ ^/^t- 
Hence G{T^y) ~ hG/R^ = {Lp/Rs)'^R~'^, which is much less than the background 
curvature provided Rg ^ Lp, i.e. provided the black hole is large compared to the 
Planck length. 

Although a tiny effect for astrophysical black holes, the Hawking process has a pro- 
found implication: black holes can "evaporate" . How long does one take to evaporate? It 
emits roughly one Hawking quantum per light crossing time, hence dM/dt ~ Th/Rs ~ 
h/Rl ~ h/G^M"^. In Planck units ;i = c = G = l, we have dM/dt ~ M'^. Integra- 
tion yields the lifetime ~ ~ {Rg/ Lp^Rs- For the 10^^ gm, 10~^^ cm black hole 
mentioned earlier we have (10^°)^10^^^cm = 10^'' cm, which is the age of the universe. 
Hence a black hole of that mass in the early universe would be explosively ending its 
life now. These are called "primordial black holes" (pbh's). None have been knowingly 
observed so far, nor their detritus, which puts limits on the present density of pbh's (see 
for example the references in^H])- However, it has been suggested [ISJ that they might 
nevertheless be the source of the highest energy (~ 3 x 10^° eV) cosmic rays. Note that 
even if pbh's were copiously produced in the early universe, their initial number density 
could easily have been inflated away if they formed before inflation. 

2.3 Black hole entropy 

How about the black hole entropy? The thermodynamic relation dM = T^dS = 
{h/8'n'GM)dS (for a black hole with no angular momentum or charge) implies Sbh = 
AirGM'^/h = Ah/4:L1,, where Ah = AirRg is the black hole horizon area. What is this 
huge entropy? What microstates does it count? 

Whatever the microstates may be, Sbh is the lower bound for the entropy emitted 
to the outside world as the black hole evaporates, i.e. the minimal missing information 
brought about by the presence of the black hole. By sending energy into a black hole 
one could make it last forever, emitting an arbitrarily large amount of entropy, so it 
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only makes sense to talk about the minimal entropy. This lower bound is attained when 
the black hole evaporates reversibly into a thermal bath at a temperature infinitesimally 
below Tfj- (Evaporation into vacuum is irreversible, so the total entropy of the outside 
increases even morejin], >S'emitted ~ {4:/3)Sbh-) This amounts to a huge entropy increase. 
The only way out of this conclusion would be if the semi-classical analysis breaks down 
for some reason... which is suspected by many, not including me. 

The so-called "information paradox" refers to the loss of information associated with 
this entropy increase when a black hole evaporates, as well as to the loss of any other 
information that falls into the black hole. I consider it no paradox at all, but many 
think it is a problem to be avoided at all costs. I think this viewpoint results from 
missing the strongly non-perturbative role of quantum gravity in evolving the spacetime 
into and beyond the classically singular region, whether by generating baby universes or 
otherwise (see section 15.3. 121 for more discussion of this point). Unfortunately it appears 
unlikely that semi-classical qg can "prove" that information is or is not lost, though 
people have tried hard. For it not to be lost there would have to be subtle effects not 
captured by semi-classical qg, yet in a regime where semi-classical description "should" 
be the whole story at leading order. 

3 Harmonic oscillator 

One can study a lot of interesting issues with free fields in curved spacetime, and I will 
restrict entirely to that case except for a brief discussion. Free fields in curved spacetime 
are similar to collections of harmonic oscillators with time-dependent frequencies. I will 
therefore begin by developing the properties of a one-dimensional quantum harmonic 
oscillator in a formalism parallel to that used in quantum field theory. 

The action for a particle of mass m moving in a time dependent potential V{x,t) in 
one dimension takes the form 



from which follows the equation of motion mx = —dxV{x,t). Canonical quantization 
proceeds by (i) defining the momentum conjugate to x, p = dL/dx = mx, (ii) replacing 
X and p by operators x and p, and (iii) imposing the canonical commutation relations 
[x,p] = ih. The operators are represented as hermitian linear operators on a Hilbert 
space, the hermiticity ensuring that their spectrum is real as befits a quantity whose 
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classical correspondent is a real number.^ In the Schrodinger picture the state is time- 
dependent, and the operators are time-independent. In the position representation for 
example the momentum operator is given hj p = —ihdx- In the Heisenberg picture the 
state is time-independent while the operators are time-dependent. The commutation 
relation then should hold at each time, but this is still really only one commutation 
relation since the the equation of motion implies that if it holds at one initial time it 
will hold at all times. In terms of the position and velocity, the commutation relation(s) 
in the Heisenberg picture take the form 

[x{t),x{t)] = ih/m. (3.2) 

Here and from here on the hats distinguishing numbers from operators are dropped. 

Specializing now to a harmonic oscillator potential V{x, t) = |mco'^(t)x^ the equation 
of motion takes the form 

x + uj^{t)x = 0. (3.3) 

Consider now any operator solution x(t) to this equation. Since the equation is second 
order the solution is determined by the two hermitian operators x{0) and x(0), and since 
the equation is linear the solution is linear in these operators. It is convenient to trade 
the pair x{0) and x{0) for a single time-independent non- hermitian operator a, in terms 
of which the solution is written as 

xit)=fit)a + ma\ (3.4) 

where f{t) is a complex function satisfying the classical equation of motion, 

f + cu'it)f = 0, (3.5) 

/ is the complex conjugate of /, and is the hermitian conjugate of a. The commutation 
relations ()3.2|) take the form 

(/,/)[a,at] = l, (3.6) 
where the bracket notation is defined by 

(/, g) = (zm/h) {fdtg - {dj)g) . (3.7) 

^In a more abstract algebraic approach, one does not require at this stage a representation but 
rather requires that the quantum variables are elements of an algebra equipped with a star operation 
satisfying certain axioms. For quantum mechanics the algebraic approach is no different from the 
concrete representation approach, however for quantum fields the more general algebraic approach 
turns out to be necessary to have sufficient generality. In these lectures I will ignore this distinction. 
For an introduction to the algebraic approach in the context of quantum fields in curved spacetime see 
|S] . For a comprehensive treatment of the algebraic approach to quantum field theory see [201 • 
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If the functions / and g are solutions to the harmonic oscillator equation ()3.5|) . then 
the bracket ()3.7j) is independent of the time t at which the right hand side is evaluated, 
which is consistent with the assumed time independence of a. 

Let us now assume that the solution / is chosen so that the real number (/, /) is 
positive. Then by rescaling / we can arrange to have 

(/,/) = !• (3.8) 

In this case the commutation relation ()3.6|) becomes 

[a,a^] = l, (3.9) 

the standard relation for the harmonic oscillator raising and lowering operators. Using 
the bracket with the operator x we can pluck out the raising and lowering operators 
from the position operator, 

a = (/,x), at = -(/,x). (3.10) 

Since both / and x satisfy the equation of motion, the brackets in ()3.10p are time 
independent as they must be. 

A Hilbert space representation of the operators can be built by introducing a state 
|0) defined to be normalized and satisfying a|0) = 0. For each n, the state \n) = 
(l/^A^^)(a^)"|0) is a normalized eigenstate of the number operator = a^a with eigen- 
value n. The span of all these states defines a Hilbert space of "excitations" above the 
state |0). 

So far the solution f{t) is arbitrary, except for the normalization condition (j3.8j) . 
A change in f{t) could be accompanied by a change in a that keeps the solution x{t) 
unchanged. In the special case of a constant frequency uit) = u however, the energy is 
conserved, and a special choice of f(t) is selected if we require that the state |0) be the 
ground state of the Hamiltonian. Let us see how this comes about. 

For a general / we have 

H = imi;2 + Imoo'^x'^ (3.11) 
'(/' + uj^P)aa + {p + uo^pfa^a^ + (|/|2 + cu^HM + •(3-12) 



Thus 



H\Q) = lm{p + u'pya^a^O) + {\ff + ooVl'm, (3.13) 
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where the commutation relation ()3.9p was used in the last term. If |0) is to be an 
eigenstate of if, the first term must vanish, which requires 



/ = ±^ujf. (3.14) 

For such an / the norm is 

(/./> = T^l/P. (3.15) 

hence the positivity of the normalization condition ()3.8p selects from ()3.14|) the minus 
sign. This yields what is called the normalized positive frequency solution to the equation 
of motion, defined by 

up to an arbitrary constant phase factor. 

With / given by ()3.16p the Hamiltonian ()3.12|) becomes 

H = lhuj{aa^ + a^a) (3.17) 
= huj{N + ^), (3.18) 

where the commutation relation (j3.9|) was used in the last step. The spectrum of the 
number operator is the non-negative integers, hence the minimum energy state is the 
one with = 0, and "zero-point energy" hjjj/2. This is just the state |0) annihilated by 
a as defined above. If any function other than ()3.16|) is chosen to expand the position 
operator as in ()3.4j) . the state annihilated by a is not the ground state of the oscillator. 

Note that although the mean value of the position is zero in the ground state, the 
mean of its square is 

(01x^10) = h/2muj. (3.19) 
This characterizes the "zero-point fluctuations" of the position in the ground state. 



4 Quantum scalar field in curved spacetime 

Much of interest can be done with a scalar field, so it suffices for an introduction. 
The basic concepts and methods extend straightforwardly to tensor and spinor fields. 
To being with let's take a spacetime of arbitrary dimension D, with a metric g^v of 
signature (H — ■■■—). The action for the scalar field is 

S= j d^x [g^^'d.^d,^ - {rr? + iR)^-") , (4.1) 
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for which the equation of motion is 

{D + m^ + ^R)y, = 0, □ = Igl-'^Xlgl'^'g^'d,. (4.2) 

(With h exphcit, the mass m should be replaced by m/h, however we'll leave h implicit 
here.) The case where the coupling ^ to the Ricci scalar R vanishes is referred to as 
"minimal coupling", and that equation is called the Klein-Gordon (KG) equation. If 
also the mass m vanishes it is called the "massless, minimally coupled scalar" . Another 
special case of interest is "conformal coupling" with m = and C, = {D — 2)/4{D — 1). 

4.1 Conformal coupling 

Let me pause briefly to explain the meaning of conformal coupling since it comes up 
often in discussions of quantum fields in curved spacetime, primarily either because 
Robertson- Walker metrics are conformally fiat or because all two-dimensional metrics 
are conformally fiat. Consider making a position dependent conformal transformation 
of the metric: 

5^^^ = n'^{x)gf,^, (4.3) 

which induces the changes 

= n-\x)g^'', \g\'^' = n^{x)\g\'/', Igl'/'r" = ^^^-'(a;)|^7^/V^ (4.4) 

R = g^^R^^ = ^-'^(r - 2{D - l)n\xin -{D-1){D~ 2)^"^(lnr]),„(lnn),^). (4.5) 

In D = 2 dimensions, the action is simply invariant in the massless, minimally coupled 
case without any change of the scalar field: S[ip,g] = S[(p,])\. In any other dimension, 
the kinetic term is invariant under a conformal transformation with a constant Q if we ac- 
company the metric change ()4.3p with a change of the scalar field, (p = QP'~^'^/'^Lp. (This 
scaling relation corresponds to the fact that the scalar field has dimension [length] 
since the action must be dimensionless after factoring out an overall h.) For a non- 
constant Vt the derivatives in the kinetic term ruin the invariance in general. However 
it can be shown that the action is invariant (up to a boundary term) if the coupling 
constant ^ is chosen to have the special value given in the previous paragraph, i.e. 
S[ip, g] = S[ip,'g\. In D = 4 dimensions that value is ,^ = 1/6. 
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4.2 Canonical quantization 

To canonically quantize we first pass to the Hamiltonian description. Separating out a 
time coordinate x'^, = {x^, x*), we can write the action as 

S = j dx^ L, L = j d^~\ C. (4.6) 

The canonical momentum at a time x^ is given by 



5\doif{x) 

Here x labels a point on a surface of constant the x^ argument of is suppressed, 
is the unit normal to the surface, and h is the determinant of the induced spatial 
metric hij. To quantize, the field and its conjugate momentum are now promoted 
to hermitian operators^ and required to satisfy the canonical commutation relation, 

Mx),'n{y)]=iM''-\x,y) (4.8) 

It is worth noting that, being a variational derivative, the conjugate momentum is a 
density of weight one, and hence the Dirac delta function on the right hand side of 
()4.8p is a density of weight one in the second argument. It is defined by the property 
/ d^~^y {x, y) f (y) = f{x) for any scalar function /, without the use of a metric 
volume element. 

In analogy with the bracket ()3.7p defined for the case of the harmonic oscillator, one 
can form a conserved bracket from two complex solutions to the scalar wave equation 

{f,g) = f{f,g) = {t/hM/'g^''(jd,g-{dJ)gy (4.9) 

This bracket is sometimes called the Klein- Gordon inner product, and (/, /) the Klein- 
Gordon norm of /. The current density j^{f,g) is divergenceless (5^j^ = 0) when the 
functions / and g satisfy the KG equation ()4.2j) , hence the value of the integral in ()4.9p is 
independent of the spacelike surface S over which it is evaluated, provided the functions 
vanish at spatial infinity. The KG inner product satisfies the relations 



{f,9) = -{f,9) = {9J), (/,/) = (4.10) 

Note that it is not positive definite. 
^See footnote^ 
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4.3 Hilbert space 



At this point it is common to expand the field operator in modes and to associate 
annihilation and creation operators with modes, in close analogy with the harmonic 
oscillator ()3.4|1 . however instead I will begin with individual wave packet solutions. My 
reason is that in some situations there is no particularly natural set of modes, and none is 
needed to make physical predictions from the theory. (An illustration of this statement 
will be given in our treatment of the Hawking effect.) A mode decomposition is a basis 
in the space of solutions, and has no fundamental status. 

In analogy with the harmonic oscillator case we define the annihilation op- 

erator associated with a complex classical solution / by the bracket of / with the field 
operator ip: 

a(/) = (/,^) (4.11) 

Since both / and 93 satisfy the wave equation, a(/) is well-defined, independent of the 
surface on which the bracket integral is evaluated. It follows from the above definition 
and the hermiticity of that the hermitian conjugate of a(/) is given by 

= -«(/)• (4.12) 

The canonical commutation relation ()4.8p together with the definition of the momentum 
PTfjl imply that 

[a(/),at(^)] = (/,^). (4.13) 

The converse is also true, in the sense that if ()4.13p holds for all solutions / and (7, 
then the canonical commutation relation holds. Using ()4.12j) . we immediately obtain 
the similar relations 

W),a{g)\ = -{f,g), [a\f),a\g)] = -{f,g) (4.14) 

If / is a positive norm solution with unit norm (/, /) = 1, then a{f) and a^(/) satisfy 
the usual commutation relation for the raising and lowering operators for a harmonic 
oscillator, [a(/), a"''(/)] = 1. Suppose now that |\I') is a normalized quantum state satis- 
fying a(/)|\I') = 0. Of course this condition does not specify the state, but rather only 
one aspect of the state. Nevertheless, for each n, the state |n, \E') = (l/v^)(a^(/))"|^E') 
is a normalized eigenstate of the number operator A^(/) = a^{f)a{f) with eigenvalue n. 
The span of all these states defines a Fock space of /-wavepacket "n-particle excitations" 
above the state l^'). 

If we want to construct the full Hilbert space of the field theory, how can we proceed? 
We should find a decomposition of the space of complex solutions to the wave equation 
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S into a direct sum of a positive norm subspace Sp and its complex conjugate Sp, such 
that all brackets between solutions from the two subspaces vanish. That is, we must 
find a direct sum decomposition 

S = Sp®S:p (4.15) 

such that 

(/,/) > \/feSp (4.16) 

{f,g) = yf,geSp. (4.17) 

The first condition implies that each / in Sp can be scaled to define its own harmonic 
oscillator sub-albegra as in the previous paragraph. The second condition implies, ac- 
cording to ()4.14j) . that the annihilators and creators for / and g in the subspace Sp 
commute amongst themselves: [a{f),a{g)] = = [a'^{f),a^{g)]. 

Given such a decompostion a total Hilbert space for the field theory can be de- 
fined as the space of finite norm sums of possibly infinitely many states of the form 
a^ifi) ■ ■ ■a^(/„)|0), where |0) is a state such that a(/)|0) = for all f in Sp, and all 
/i, . . . , /„ are in Sp. The state |0) is called a Fock vacuum. It depends on the decom- 
position ()4.15|) . and in general is not the ground state (which is not even well defined 
unless the background metric is globally static). The representation of the field operator 
on this Fock space is hermitian and satisfies the canonical commutation relations. 

4.4 Flat spacetime 

Now let's apply the above generalities to the case of a massive scalar field in flat space- 
time. In this setting a natural decomposition of the space of solutions is defined by 
positive and negative frequency with respect to a Minkowski time translation, and the 
corresponding Fock vacuum is the ground state. I summarize briefly since this is stan- 
dard fiat spacetime quantum field theory. 

Because of the infinite volume of space, plane wave solutions are of course not nor- 
malizable. To keep the physics straight and the language simple it is helpful to introduce 
periodic boundary conditions, so that space becomes a large three-dimensional torus with 
circumferences L and volume V = L^. The allowed wave vectors are then k = (27r/L)n, 
where the components of the vector n are integers. In the end we can always take the 
limit L ^ oo to obtain results for local quantities that are insensitive to this formal 
compactification. 



15 



A complete set of solutions ("modes") to the classical wave equation ()4.2|) with □ 
the fiat space d'Alembertian and = is given by 



where 

uj{k) = Vk2 + m2, (4.19) 

together with the solutions obtained by replacing the positive frequency u;(k) by its 
negative, —uj(k). The brackets between these solutions satisfy 

(/k,/i) = 4,1 (4.20) 
(7k, Ti) = -Vi (4.21) 
(/k,7i) = 0, (4.22) 

so they provide an orthogonal decomposition of the solution space into positive norm 
solutions and their conjugates as in ()4.17p . with Sp the space spanned by the positive 
frequency modes f^. As described in the previous subsection this provides a Fock space 
representation. 

If we define the annihilation operator associated to /k by 

ak=(/k,<^), (4.23) 
then the field operator has the expansion 

^ = Y1 (/kak+7kal)- (4.24) 

k 

Since the individual solutions /k have positive frequency, and the Hamiltonian is a sum 
over the contributions from each k value, our previous discussion of the single oscillator 
shows that the vacuum state defined by 

ak|0) = (4.25) 

for all k is in fact the ground state of the Hamiltonian. The states 

4|0) (4.26) 

have momentum ^k and energy /ia;(k), and are interpreted as single particle states. 
States of the form aj^^ ■ ■ ■ a|^^|0) are interpreted as n-particle states. 
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Note that although the field Fourier component y^k = /k Ok + /_k '^-k ^^^^ 
mean in the vacuum state, like the harmonic oscillator position it undergoes "zero-point 
fiuctuations" characterized by 



which is entirely analogous to the oscillator result ()3.19|) . 

4.5 Curved spacetime, "particles", and stress tensor 

In a general curved spacetime setting there is no analog of the preferred Minkowski 
vacuum and definition of particle states. However, it is clear that we can import these 
notions locally in an approximate sense if the wavevector and frequency are high enough 
compared to the inverse radius of curvature. Slightly more precisely, we can expand the 
metric in Riemann normal coordinates about any point Xq: 



If and i-i;^(k) are much larger than any component of the Riemann tensor R^i,ai3{xo) 
in this coordinate system then it is clear that the fiat space interpretation of the cor- 
responding part of Fock space will hold to a good approximation, and in particular a 
particle detector will respond to the Fock states as it would in fiat spacetime. This 
notion is useful at high enough wave vectors locally in any spacetime, and for essentially 
all wave vectors asymptotically in spacetimes that are asymptotically fiat in the past or 
the future or both. 

More generally, however, the notion of a "particle" is ambiguous in curved space- 
time, and one should use field observables to characterize states. One such observable 
determines how a "particle detector" coupled to the field would respond were it following 
some particular worldline in spacetime and the coupling were adiabatically turned on 
and off at prescribed times. For example the transition probability of a point monopole 
detector is determined in lowest order perturbation theory by the two-point function 
{'^\ip{x)Lp{x')\'^) evaluated along the worldhnePTl ITT] of the detector. (For a careful 
discussion of the regular ization required in the case of a point detector see |^.) Al- 
ternatively this quantity — along with the higher order correlation functions — is itself a 
probe of the state of the field. 

Another example of a field observable is the expectation value of the stress energy 
tensor, which is the source term in the semi-classical Einstein equation 



{o\ipW^\o) = |/-k 



h 



(4.27) 



2Vuj{ky 



9^,u{x) = r]^iu + -R^^ua|3{xo){x - XoT{x - XqY + 0{{x - x^f). 





(4.29) 
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This quantity is infinite because it contains the product of field operators at the same 
point. Physically, the infinity is due to the fluctuations of the infinitely many ultraviolet 
field modes. For example, the leading order divergence of the energy density can be 
attributed to the zero-point energy of the field fluctuations, but there are subleading 
divergences as well. We have no time here to properly go into this subject, but rather 
settle for a few brief comments. 

One way to make sense of the expectation value is via the difference between its 
values in two different states. This difference is well-defined (with a suitable regulator) 
and finite for any two states sharing the same singular short distance behavior. The 
result depends of course upon the comparison state. Since the divergence is associated 
with the very short wavelength modes, it might seem that to uniquely define the ex- 
pectation value at a point x it should suffice to just subtract the infinities for a state 
defined as the vacuum in the local flat spacetime approximation at x. This subtraction 
is ambiguous however, even after making it as local as possible and ensuring local con- 
servation of energy {y^^{T^y) = 0). It defines the expectation value only up to a tensor 
H^j^y constructed locally from the background metric with four or fewer derivatives and 
satisfying the identity V^H^^y = 0. The general such tensor is the variation with respect 
to the metric of the invariant functional 

J rf^x (co + ciR + C2R^ + c^R^'R^u + c^R^"'"' R^.upa) . (4.30) 

(In four spacetime dimensions the last term can be rewritten as a combination of the 
first two and a total divergence.) Thus H^y is a combination of (7^,^, the Einstein tensor 
G^jy, and curvature squared terms. In effect, the ambiguity —H^^ is added to the metric 
side of the semi-classical field equation, where it renormalizes the cosmological constant 
and Newton's constant, and introduces curvature squared terms. 

A different approach is to define the expectation value of the stress tensor via the 
metric variation of the renormalized effective action, which posesses ambiguities of the 
same form as ()4.3U|) . Hence the two approaches agree. 

4.6 Remarks 

4.6.1 Continuum normalization of modes 

Instead of the "box normalization" used above we could "normalize" the solutions /k 
(prTR|l with the factor V^^/^ replaced by (27r)-3/2. Then the Kronecker 5's in (jODllOII) 
would be replaced by Dirac 5-functions and the discrete sum over momenta in ()4.24|) 
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would be an integral over k. In this case the annihilation and creation operators would 
satisfy [flk, ctj] = ^'^(k, I). 

4.6.2 Massless minimally coupled zero mode 

The massless minimally coupled case m = 0, ^ = has a peculiar feature. The spatially 
constant function f{t, x) = Co + Cit is a solution to the wave equation that is not included 
among the positive frequency solutions /k or their conjugates. This "zero mode" must 
be quantized as well, but it behaves like a free particle rather than like a harmonic 
oscillator. In particular, the state of lowest energy would have vanishing conjugate field 
momentum and hence would be described by a Schrodinger wave function ip{(po) that 
is totally delocalized in the field amplitude ipo. Such a wave function would be non- 
normalizable as a quantum state, just as any momentum eigenstate of a non-relativistic 
particle is non-normalizable. Any normalized state would be described by a Schrodinger 
wavepacket that would spread in (po like a free particle spreads in position space, and 
would have an expectation value ((/Pq) growing linearly in time. This suggests that no 
time-independent state exists. That is indeed true in the case where the spatial directions 
are compactified on a torus for example, so that the modes are discrete and the zero 
mode carries as much weight as any other mode. In non-compact space one must look 
more closely. It turns out that in 1 + 1 dimensions the zero mode continues to preclude a 
time independent state (see e.g. although the connection to the behavior of the zero 
mode is not made there), however in higher dimensions it does not, presumably because 
of the extra factors of k in the measure k^~'^dk. A version of the same issue arises 
in deSitter space, where no deSitter invariant state exists for the massless, minimally 
coupled field That is true in higher dimensions as well, however, which may be 

related to the fact that the spatial sections of deSitter space are compact. 

5 Particle creation 

We turn now to the subject of particle creation in curved spacetime (which would more 
appropriately be called "field excitation", but we use the standard term). The main 
applications are to the case of expanding cosmological spacetimes and to the Hawking 
effect for black holes. To begin with I will discuss the analogous effect for a single har- 
monic oscillator, which already contains the essential elements of the more complicated 
cases. The results will then be carried over to the cosmological setting. The following 
section then takes up the subject of the Hawking effect. 
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5.1 Parametric excitation of a harmonic oscillator 



A quantum field in a time-dependent background spacetime can be modeled in a simple 
way by a harmonic oscillator whose frequency uj{t) is a given function of time. The 
equation of motion is then Q, 

x + uj^{t)x = 0. (5.1) 

We consider the situation where the frequency is asymptotically constant, approaching 
Uin in the past and Uout in the future. The question to be answered is this: if the 
oscillator starts out in the ground state |0j„) appropriate to as t — — oo, what is the 
state as t ^ +00? More precisely, in the Heisenberg picture the state does not evolve, 
so what we are really asking is how is the state |0j„) expressed as a Fock state in the 
out-Hilbert space appropriate to tUo^t? It is evidently not the same as the ground state 
\Oout) of the Hamiltonian in the asymptotic future. 

To answer this question we need only relate the annihilation and creation operators 
associated with the in and out normalized positive frequency modes f in (t), which are 

out 

solutions to (jS.lj) with the asymptotic behavior 



fin{t) > eXp{-%UJ in t) . (5.2) 

out \ ZniLU in out 

V out 

Since the equation of motion is second order in time derivatives it admits a two-parameter 
family of solutions, hence there must exist complex constants a and /3 such that 

fout = a fin + Pfin- (5.3) 

The normalization condition {fout, fout) = 1 implies that 

|«|2-|/?|2 = l. (5.4) 
The out annihilation operator is given (see ()3.10|) ) by 

(^out — ifoutj^) (5-5) 
= {a fin + (3 fin, x) (5.6) 

= aain - (3a\^. (5.7) 

This sort of linear relation between two sets of annihilation and creation operators (or the 
corresponding solutions) is called a Bogoliuhov transformation, and the coefficients are 
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the Bogoliuhov coefficients. The mean value of the out number operator Nout = ciout^'out 
is nonzero in the state |Ojn): 

{Qin\Nout%n) = \I3\\ (5.8) 

In this sense the time dependence of uj{t) excites the oscillator, and characterizes 
the excitation number. 

To get a feel for the Bogoliubov coefficient let us consider two extreme cases, 
adiabatic and sudden. 

5.1.1 Adiabatic transitions and ground state 

The adiabatic case corresponds to a situation in which the frequency is changing very 
slowly compared to the period of oscillation, 

-<cj. 5.9 

In this case there is almost no excitation, so \(3\ ^ 1. Generically in the adiabatic case 
the Bogoliubov coefficient is exponentially small, [3 ~ exp(— cuoT"), where ujq is a typical 
frequency and T characterizes the time scale for the variations in the frequency. 

In the Schrodinger picture, one can say that during an adiabatic change of uj{t) 
the state continually adjusts to remain close to the instantaneous adiabatic ground state. 
This state at time to is the one annihilated by the lowering operator a^fto) defined by the 
solution (t) satisfying the initial conditions at to corresponding to the "instantaneous 
positive frequency solution", 

^(to) = ^h/2muj{to) (5.10) 
hM = -iuj{t^)ftM. (5.11) 

The instantaneous adiabatic ground state is sometimes called the "lowest order adiabatic 
ground state at to" . One can also consider higher order adiabatic ground states as follows 
(see e.g. ^^). A function of the WKB form {h/2mW{t)f''^ ey.Y>{-i !^W{t')dt') is a 
normalized solution to ()5.H1 provided W{t) satisfies a certain second order differential 
equation. That equation can be solved iteratively, yielding an expansion W{t) = uj{t) + 
■ ■ -, where the subsequent terms involve time derivatives of uj{t). The lowest order 
adiabatic ground state at to is defined using the solution whose initial conditions ()5.11|1 
match the lowest order truncation of the expansion for W{t). A higher order adiabatic 
ground state is similarly defined using a higher order truncation. 
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5.1.2 Sudden transitions 



The opposite extreme is tlie sudden one, in which u changes instantaneously from Um 
to Uout at some time to- We can then find the Bogohubov coefficients using ()5.3|) and its 
first derivative at to- For to = the result is 



2 V V ^out V 



a = Ua/^+a/^) (5-12) 



P = U^F^-^F^]■ (5.13) 



2 V V ^out V 

(For to 7^ there are extra phase factors in these solutions.) Interestingly, the amount 
of excitation is precisely the same if the roles of and Uout are interchanged. For an 
example, consider the case where Uout = '^^mi for which a = 5/4 and (3 = —3/4. In this 
case the expectation value ()5.8|1 of the out number operator is 9/16, so there is about 
"half an excitation" . 

5.1.3 Relation between in and out ground states &; the squeeze operator 

The expectation value of Nout is only one number characterizing the relation between 
\Oin) and the out states. We shall now determine the complete relation 

\Oin) =^Cn\n)out (5.14) 
n 

where the states \n)out are eigenstates of Nout and the c„ are constants. 

One can find am in terms of aout and a^^^ by combining ()5.3|1 with its complex 
conjugate to solve for fin in terms of fout and fout- In analogy with ()5.7|) one then finds 

Qin = aaout + f^ciiuf (5.15) 
Thus the defining condition aj„|Om) = implies 

aout\^in) = -- a^out\^in) ■ (5.16) 

a 

A transparent way to solve this is to note that the commutation relation [aout^ c^Lj = 1 
suggests the formal analogy aout = d/da\^^. This casts (j5.16|) as a first order ordinary 
differential equation, with solution 



|0m) = ATexp 



2a 



^out^out 



|0„„,) (5.17) 
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Note that the state \Oout) on which the exponential operator acts is annihilated by aout, 
so there is no extra term from aou^-dependence. 

The state |0m) thus contains only even numbered excitations when expressed in terms 
of the out number eigenstates. The normalization constant A/" is given by 



2n! 



2a 



2n 



(5.19) 



For large n the summand approaches l/^/ap" = + where the relation 

()5.4p was used in the last step. The sum therefore converges, and can be evaluated^ to 
yield 



lATI 



\P/a\ 



1/4 



\a\ 



-1/2 



(5.20) 



An alternate way of describing |0j„) in the out Hilbert space is via the squeeze operator 



S = exp 



— a'a' — —aa 



.2 2 

Since its exponent is anti-hermitian, S is unitary. Conjugating a by S* yields 

z + 



S^aS = cosh \z\ a + sinh 



z\ 



a 



\z\ 



(5.21) 



(5.22) 



This has the form of the Bogoliubov transformation ()5.15|1 with a = cosh \ z\ and (3 
sinh I2; 1(^/1^1). With aout in place of a in 5* this gives 



Oiin — S^aQutS. 

The condition ai„|Om) = thus implies aoutS\Oin) = 0, so evidently 

|0..) = S^Oout) 



(5.23) 



(5.24) 



up to a constant phase factor. That is, the in and out ground states are related by the 
action of the squeeze operator S. Since S is unitary, the right hand side of ()5.24|) is 
manifestly normalized. 



■^Alex Maloney pointed out that it can be evaluated by expressing the binomial coefhcient 2n!/(n!)^ 
as the contour integral f{dz/2'!riz){z+ and interchanging the order of the sum with the integral. 
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5.2 Cosmological particle creation 

We now apply the ideas just developed to a free scalar quantum field satisfying the KG 
equation ()4.2j) in a homogeneous isotropic spacetime. The case when the spatial sections 
are flat is slightly simpler, and it is already quite apphcable, hence we restrict to that 
case here. 

The spatially flat Robertson- Walker (RW) line element takes the form: 

ds'^ = dt^ - a^{t)dx'dx' = a^{i]){dif - dx'dx') (5.25) 

It is conformally fiat, as are all RW metrics. The coordinate f] = J dt/a{t) is called the 
conformal time, to distinguish it from the proper time t of the isotropic observers. The 
d'Alembertian □ for this metric is given by 

n = d^t +{?>d/a)dt-a~^dl, (5.26) 

where the dot stands for d/dt. The spatial translation symmetry allows the spatial 
dependence to be separated from the time dependence. A field 

Mfc(x,t) = Ck(t)e''^'' (5.27) 

satisfies the field equation ()4.2|) provided Ck(^) satisfies an equation similar to that of 
a damped harmonic oscillator but with time-dependent damping coefficient 3a/a and 
time-dependent frequency a~^/c^ -|-m^ + ^R- It is worth emphasizing that in spite of the 
"damping", the field equation is Hamiltonian, and the Klein-Gordon norm ()4.9|) of any 
solution is conserved in the evolution. 

The field equation can be put into the form of an undamped oscillator with time- 
dependent frequency by using the conformal time 7] instead of t and factoring out an 
appropriate power of the conformal factor a?{ri): 

Ck = a-iXk. (5.28) 

The function satisfies the field equation if and only if 

x;; + a;2(ry)Xk = 0, (5.29) 

where the prime stands for d/dr] and 

uj\r]) = e + m^a^ - (1 - G^W'/a). (5.30) 
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(In the special case of conformal coupling m = and ^ = 1/6, this becomes the time- 
independent harmonic oscillator, so that case is just like flat spacetime. All effects of 
the curvature are then incorporated by the prefactor a{r])~^ in (j5.28|) .) 

The Mk are orthogonal in the Klein-Gordon inner product ()4.9|) . and they are nor- 
malized^ provided the Xk have unit norm: 

(uk,Mi) = 5k,i iiV/h){xkXk - XkXk) = 1- (5.31) 

Note that the relevant norm for the Xk differs from that for the harmonic oscillator 
()3.7p only by the replacement m ^ V, where V is the a;*-coordinate volume of the 
constant t surfaces. We also have (uk, ui) = for all k, 1, hence these modes provide 
an orthogonal postive/negative norm decomposition of the space of complex solutions. 
As discussed in section 14.31 this yields a corresponding Fock space representation for 
the field operators. The field operator can be expanded in terms of the corresponding 
annihilation and creation operators: 

V?(x, t) =^ (^Mk(x, t) ak + Mk(x, t) al^ . (5.32) 

k 

Consider now the special case where there is no time dependence in the past and 
future, a{ri) constant. The in and out "vacua" are the ground states of the Hamilto- 
nian at early and late times, and are the states annihilated by the Ok associated with the 
^in,out QQnstructed with early and late time positive frequency modes xjf as explained 
in section 

X^\V) expi-tu; v). (5.33) 

\/ ZV U in out 

V out 

The Bogoliubov transformation now takes the form 

= E + ^kk'wic") • (5.34) 

k' 

Matching the coefficients of exp(ik ■ x), we see that 

ttkk' = ttfc^k.k', /?kk' = ^fc^k-k', (5.35) 

i.e. the Bogoliubov coefficients mix only modes of wave vectors k and — k, and they 
depend only upon the magnitude of the wavevector on account of rotational symmetry 



''The two inverse factors of a coming from (|5.28|) are cancelled by the in the volume element and 
the in the relation d/dt — a^^d/drj 
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(eqn. ()5.29|) for does not depend on the direction of k). The normahzation condition 
on imphes 

- = 1. (5.36) 

As in the harmonic oscillator example ()5.8p . if the state is the m- vacuum, then the 
expected excitation level of the k out-mode, i.e. the average number of particles in that 
mode, is given by 

(0.„|iVr*|0,„) = (5.37) 

To convert this statement into one about particle density, we sum over k and divide by 
the physical spatial volume Vphys = a^V, which yields the number density of particles. 
Alternatively one can work with the continuum normalized modes. The relation between 
the discrete and continuous sums is 

77— Y ' — ' I ^^'k. (5.38) 

the number density of out-particles is thus 

The mean particle number characterizes only certain aspects of the state. As in the 
oscillator example, a full description of the in-vacuum in the out-Fock space is obtained 
from the Bogoliubov relation between the corresponding annihilation and creation op- 
erators. From (j5.34j) we can solve for m*" and thence find 

< = «fe«r* + A«V"*, (5-40) 



whence 



ar*|0.„) = -^aV"*|0,„), (5.41) 

Oik 



which can be solved to find 



l^k \ ^toMt^t out 

2ak 



\Oout) (5.42) 



k 

where the J\f{^ are normalization factors. This solution is similar to the corresponding 
expression ()5.18p for the harmonic oscillator and it can be found by a similar method. 
Although the operators aj^°"* and alj^°"* are distinct, each product appears twice in 
the sum, once for k and once for — k, hence the factor of 2 in the denominator of the 
exponent is required. Using the two-mode analog of the squeeze operator ()5.21|1 the 
state ()5.42j) can also be written in a manifestly normalized fashion analogous to ()5.24j) . 
It is sometimes called a squeezed vacuum. 
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5.3 Remarks 



5.3.1 Momentum correlations in the squeezed vacuum 

The state ()5.42j) can be expressed as a sum of terms each of which has equal numbers of 
k and — k excitations. These degrees of freedom are thus entangled in the state, in such 
a way as to ensure zero total momentum. This is required by translation invariance of 
the states |0j„) and \Oout), since momentum is the generator of space translations. 

5.3.2 Normalization of the squeezed vacuum 

The norm sum for the part of the state ()5.42|1 involving k and — k is a standard geometric 
series, which evaluates to la^p using ()5.36|) . Hence to normalize the state one can set 
A/k = for all k, including k = as in ()5.2()|1 . The overall normalization factor 

is a product of infinitely many numbers less than unity. Unless those numbers converge 
rapidly enough to unity the state is not normalizable. The condition for normalizability 
is easily seen to be J2k l/^^P < according to ()5.H7|1 the average total number 

of excitations must be finite. If it is not, the state \Oin) does not lie in the Fock space 
built on the the state \Oout)- Note that, although formally unitary, the squeeze operator 
does not act unitarily on the out Fock space if the corresponding state is not in fact 
normalizable. 

5.3.3 Energy density 

If the scale factor a changes by over a time interval At, then for a massless field dimen- 
sional analysis indicates that the in vacuum has a resulting energy density p ~ /i(Ar)~^ 
after the change. To see how the formalism produces this, according to ()5.H9|1 we have 
p = (27ra)~^ J cPk\Pk\'^{hju}/a). The Bogoliubov coefficent Pk is of order unity around 
kc/a ~ 1/Ar and decays exponentially above that. The integral is dominated by the 
upper limit and hence yields the above mentioned result. 

5.3.4 Adiabatic vacuum 

Modes with frequency much larger than a' /a see the change of the scale factor as adia- 
batic, hence they remain relatively unexcited. The state that corresponds to the instan- 
taneously defined ground state, in analogy with ()5.1H) for the single harmonic oscillator, 
is called the adiabatic vacuum at a given time. 
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5.4 de Sitter space 



The special case of de Sitter space is of interest for various reasons. The first is just 
its high degree of symmetry, which makes it a convenient arena for the study of qft in 
curved space. It is the maximally symmetric Lorentzian space with (constant) positive 
curvature. Maximal symmetry refers to the number of Killing fields, which is the same 
as for flat spacetime. The Euclidean version of de Sitter space is just the sphere. 

de Sitter (dS) space has hypersurface-orthogonal timelike Killing fields, hence is lo- 
cally static, which further simplifies matters, but not to the point of triviality. The reason 
is that all such Killing fields have Killing horizons, null surfaces to which they are tan- 
gent, and beyond which they are spacelike. Hence dS space serves as a highly symmetric 
analog of a black hole spacetime. In particular, a symmetric variant of the Hawking 
effect takes place in de Sitter space, as first noticed by Gibbons and Hawking ^H]. See 
for a recent review. 

Inflationary cosmology provides another important use of deSitter space, since dur- 
ing the period of exponential expansion the spacetime metric is well described by dS 
space. In this application the dS line element is usually written using spatially flat RW 
coordinates: 

ds^ = df - e^^'dx'dxK (5.43) 

These coordinates cover only half of the global dS space, and they do not make the 
existence of a time translation symmetry manifest. This takes the conformal form ()5.25|) 
with 7] = —H~^ exp{—Ht) and a{ri) = —1/Ht]. The range of t is (— oo, oo) while that of 
T] is (— oo, 0). 

The flat patch of de Sitter space is asymptotically static with respect to conformal 
time 7] in the past, since a' /a = —1/ t] ^ a.s rj —>■ — oo. Therefore in the asymptotic 
past the adiabatic vacuum (with respect to positive 77-frequency) defines a natural initial 
state. This is the initial state used in cosmology. In fact it happens to define a deSitter 
invariant state, also known as the Euclidean vacuum or the Bunch-Davies vacuum. 

5.4.1 Primordial perturbations from zero point fluctuations 

Observations of the Cosmic Microwave Background radiation support the notion that 
the origin of primordial perturbations lies in the quantum fluctuations of scalar and 
tensor metric modes (see for a recent review and [21] for a classic reference.) The 
scalar modes arise from (and indeed are entirely determined by, since the metric has no 
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independent scalar degree of freedom) coupling to matter.^ Let's briefly discuss how this 
works for a massless minimally coupled scalar field, which is just how these perturbations 
are described. First I'll describe the scenario in words, then add a few equations. 

Consider a field mode with a frequency high compared to the expansion rate a/ a 
during the early universe. To be specific let us assume this rate to be a constant if, 
i.e. de Sitter inflation. Such a mode was presumably in its ground state, as the prior 
expansion would have redshifted away any initial excitation. As the universe expanded 
the frequency redshifted until it became comparable to the expansion rate, at which 
point the oscillations ceased and the field amplitude approached a time-independent 
value. Just before it stopped oscillating the field had quantum zero point fluctuations 
of its amplitude, which were then preserved during the further expansion. Since the 
amplitude was frozen when the mode had the fixed proper wavenumber if, it is the 
same for all modes apart from the proper volume factor in the mode normalization 
which varies with the cosmological time of freezeout. Finally after inflation ended the 
expansion rate dropped faster than the wavenumber, hence eventually the mode could 
begin oscillating again when its wavelength became shorter than the Hubble length H~^. 
This provided the seeds for density perturbations that would then grow by gravitational 
interactions. On account of the particular wavevector dependence of the amplitude of 
the frozen spatial fluctuations, the spectrum of these perturbations turns out to be 
scale-invariant (when appropriately defined). 

More explicitly, the field equation for a massless minimally coupled field is given by 
dESni), with 

uj\7^) = e- a" /a = e- 2H^a\ (5.44) 

where the last equality holds in de Sitter space. In terms of proper frequency ujp = oj /a 
and proper wavenumber kp = k/a we have Up = kp — 2H^. The first term is the usual 
flat space one that produces oscillations, while the second term tends to oppose the 
oscillations. For proper wavenumbers much higher than H the second term is negligible. 
The field oscillates while the proper wavenumber redshifts exponentially. Eventually the 
two terms cancel, and the mode stops oscillating. This happens when 

kp = V2H. (5.45) 

As the wavenumber continues to redshift into the region kp <^ H, to a. good approxima- 
tion the amplitude satisfies the equation Xk ~ (c^"/'^)Xk = 0; which has a growing and a 

^The Cosmic Microwave Background observations supporting this account of primordial perturba- 
tions thus amount to quantum gravity observations of a hmited kind. 
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decaying solution^. The growing solution is Xk oc a, which imphes that the field mode 
Mk ()5.27|5.28|1 is constant in time. (This conclusion is also evident directly from the fact 
that the last term of the wave operator (j5.26|) vanishes as a grows.) 

The squared amplitude of the fiuctuations when frozen is, according to ()4.27p . 

(Ol^t^klO) = l^-^r ~ ^ = (5-46) 

where the proper values are used for consistent matching to the previous fiat space result. 
This gives rise to the scale invariant spectrum of density perturbations. 



6 Black hole evaporation 

The vacuum of a quantum field is unstable to particle emission in the presence of a black 
hole event horizon. This instability is called the Hawking effect. Unlike the cosmological 
particle creation discussed in the previous section, this effect is not the result of time 
dependence of the metric exciting the field oscillators. Rather it is more like pair creation 
in an external electric field[13j. (For more discussion of the role of time dependence see 
section 16.3.131 ) A general introduction to the Hawking effect was given in section YI7I\ 
The present section is devoted to a derivation and discussion related topics. 

The historical roots of the Hawking effect lie in the classical Penrose process for 
extracting energy from a rotating black hole. We first review that process and indicate 
how it led to Hawking's discovery. Then we turn to the qft analysis. 



6.1 Historical sketch 

The Kerr metric for a rotating black hole is stationary, but the asymptotic time transla- 
tion Killing vector x becomes spacelike outside the event horizon. The region where it is 
spacelike is called the ergoregion. The conserved Killing energy for a particle with four- 
momentum p is E = X ' P- Physical particles have future pointing timelike 4-momenta, 
hence E is positive provided x is also future timelike. Where x is spacelike however, 
some physical 4-momenta have negative Killing energy. 

In the Penrose process, a particle of energy i^o > is sent into the ergoregion of a 
rotating black hole where it breaks up into two pieces with Killing energies Ei and E2, 
so that Eq = El + E2. If E2 is arranged to be negative, then Ei > Eq, that is, more 
energy comes out than entered. The black hole absorbs the negative energy E2 and thus 

^The general solution is cia + C2a drj/a^, which is csrj^^ + c^if' in de Sitter space. 
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loses mass. It also loses angular momentum, hence the process in effect extracts the 
rotational energy of the black hole. 

The Penrose process is maximally efficient and reversible if the horizon area is un- 
changed. That condition is achievable in the limit that the absorbed particle enters the 
black hole on a trajectory tangent to one of the null generators of the horizon. This 
role of the horizon area in governing efficiency of energy extraction exhibits the analogy 
between area and entropy. Together with Bekenstein's information theoretic arguments, 
it gave birth to the subject of black hole thermodynamics. 

When a field scatters from a rotating black hole a version of the Penrose process 
called superradiant scattering can occur. The analogy with stimulated emission suggests 
that quantum fields should exhibit spontaneous emission from a rotating black hole. To 
calculate this emission rate from an "eternal" spinning black hole one must specify a 
condition on the state of the quantum field that determines what emerges from the past 
horizon. In order to avoid this unphysical specification Hawking considered instead a 
black hole that forms from collapse, which has no past horizon. In this case only initial 
conditions before the collapse need be specified. 

Much to his surprise. Hawking found that for any initial state, even a non-rotating 
black hole will spontaneously emit radiation. This Hawking effect is a pair creation 
process in which one member of the pair lies in the ergoregion inside the horizon and 
has negative energy, while the other member lies outside and escapes to infinity with 
positive energy. The Hawking radiation emerges in a steady fiux with a thermal spectrum 
at the temperature Th = hK,/27T. The surface gravity k had already been seen to play 
the role of temperature in the classical first law of black hole mechanics, which thus 
rather remarkably presaged the quantum Hawking effect. 

6.2 The Hawking effect 

Two different notions of "frequency" are relevant to this discussion. One is the "Killing 
frequency", which refers to time dependence with respect to the time-translation sym- 
metry of the background black hole spacetime. In the asymptotically flat region at 
infinity the Killing frequency agrees with the usual frequency defined by the Minkowski 
observers at rest with respect to the black hole. The other notion is "free-fall frequency" 
defined by an observer falling across the event horizon. Since the Killing fiow is tangent 
to the horizon. Killing frequency there is very different from free-fall frequency, and that 
distinction hes at the heart of the Hawking effect. 
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6.2.1 Average number of outgoing particles 

The question to be answered is tliis: if a black liole forms from collapse with the quantum 
field in any 'regular' state |\E'), then at late times, long after the collapse, what will be the 
average particle number and other observables for an outgoing positive Killing frequency 
wavepacket P of a quantum field far from the black hole? We address this question here 
for the case of a noninteracting scalar field and a static (non-rotating) black hole At 
the end we make some brief remarks about generalizations. Figure Q depicts the various 
ingredients in the following discussion. 



Figure 1: Spacetime diagram of black hole formed by collapsing matter. The outgoing 
wavepacket P splits into the transmitted part T and reflected part R when propagated back- 
wards in time. The two surfaces Sj^j are employed for evaluating the Klein-Gordon inner 
products between the wavepacket and the field operator. Although P, and hence R and T, 
have purely positive Killing frequency, the free-fall observer crossing T just outside the horizon 
sees both positive and negative frequency components with respect to his proper time. 

To begin with we evaluate the expectation value (^|A(P)|\E') of the number operator 



collapsing matter 
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N{P) = a\P)a{P) in the quantum state of the field. This does not fully characterize 
the state, but it will lead directly to considerations that do. 

The annihilation operator (|4.1ip corresponding to a normalized wavepacket P is 
given by 

a(P) = (P,y.)s„ (6.1) 

where the Klein-Gordon inner product ()4.9|) is evaluated on the "final" spacelike slice 
Tjf. To evaluate the expectation value of N{P) we use the field equation satisfied by (f 
to relate N{P) to an observable on an earlier slice Sj on which we know enough about 
the quantum state. Specifically, we assume there are no incoming excitations long after 
the black hole forms, and we assume that the state looks like the vacuum at very short 
distances (or high frequencies) as seen by observers falling across the event horizon. 
Hawking originally propagated the field through the time-dependent collapsing part of 
the metric and back out all the way to spatial infinity, where he assumed |\E') to be 
the incoming vacuum at very high frequencies. As pointed out by Unruhj2Tj (see also 
fl^ IHOj) the result can be obtained without propagating all the way back, but rather 
stopping on a spacelike surface Sj far to the past of S j but still after the formation of the 
black hole. This is important since propagation back out to infinity invokes arbitrarily 
high frequency modes whose behavior may not be given by the standard relativistic free 
field theory. 

If Sj lies far enough to the past of Sj, the wavepacket P propagated backwards by 
the Klein-Gordon equation breaks up into two distinct parts, 

P = R + T. (6.2) 

(See Fig. H) R is the "refiected" part that scatters from the black hole and returns 
to large radii, while T is the "transmitted" part that approaches the horizon. R has 
support only at large radii, and T has support only in a very small region just outside 
the event horizon where it oscillates very rapidly due to the backwards gravitational 
blueshift. Since both the wavepacket P and the field operator ip satisfy the Klein-Gordon 
equation, the Klein-Gordon inner product in ()6.H) can be evaluated on Sj instead of 
without changing a{P). This yields a corresponding decomposition for the annihilation 
operator, 

a{P) = a{R) + a{T). (6.3) 

Thus we have 

(^lAr(P)l^) = (*|(a^(i?) + a\T)){a{R) + a(T))|^). (6.4) 
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Now it follows from the stationarity of the black hole metric that the Killing fre- 
quencies in a solution of the KG equation are conserved. Hence the wavepackets R and 
T both have the same, purely positive Killing, frequency components as P. As R lies 
far from the black hole in the nearly flat region, this means that it has purely positive 
asymptotic Minkowski frequencies, hence the operator a{R) is a bona fide annihilation 
operator — or rather {R, R)^/"^ times the annihilation operator — for incoming excitations. 
Assuming that long after the black hole forms there are no such incoming excitations, 
we have a{R)\^) = 0. Equation ()6.4p then becomes 

{^\N{P)\^) = {^\a\T)a{T)\^). (6.5) 

If a(T) |\E') = as well, then no P-particles are emitted at all. The state with this 
property is called the "Boulware vacuum". It is the state with no positive Killing 
frequency excitations anywhere, including at the horizon. 

The Boulware vacuum does not follow from collapse however. The reason is that 
the wavepacket T does not have purely positive frequency with respect to the time of 
a free fall observer crossing the horizon, and it is this latter frequency that matches to 
the local Minkowski frequency in a neighborhood of the horizon small compared with 
the radius of curvature of the spacetime. 

More precisely, consider a free-fall observer (i.e. a timelike geodesic x{t)) with proper 
time r who falls across the horizon at r = at a point where the slice Sj meets the 
horizon (see Fig. Q). For this observer the wavepacket T has time dependence T(r) (i.e. 
T(x(r))) that vanishes for r > since the wavepacket has no support behind the horizon. 
Such a function cannot possibly have purely positive frequency components. To see why, 
recall that if a function vanishes on a continuous arc in a domain of analyticity, then it 
vanishes everywhere in that domain (since its power series vanishes identically on the 
arc and hence by analytic continuation everywhere). Any positive frequency function 

POD _ 

/i(r) = / duj e~'^^ h{uj) (6.6) 
Jo 

is analytic in the lower half r plane, since the addition of a negative imaginary part to 
r leaves the integral convergent. The positive real r axis is the limit of an arc in the 
lower half plane, hence if /i(r) were to vanish for r > it would necessarily vanish also 
for r < 0. (Conversely, a function that is analytic on the lower half-plane and does not 
blow up exponentially as |r| oo must contain only positive frequency components, 
since exp(— zcur) does blow up exponentially as |r| ^ cxd when uj is negative.) 
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The wavepacket T can be decomposed into its positive and negative frequency parts 
with respect to the free fall time r, 



(6.7) 



which yields the corresponding decomposition of the annihilation operator 



a(T) 



a{T+) + a{T-) 



(6.8) 
(6.9) 



Since T~ has negative KG norm, Eqn. ()4.12|) has been used in the last line to trade 
a(T~) for the bona fide creation operator a^{T~). The r-dependence of T consists of 
very rapid oscillations for r < 0, so the wavepackets T"*" and T~ have very high energy 
in the free-fall frame. 

A free fall observer crossing the horizon long after the black hole forms would pre- 
sumably see the ground state of the field at short distances, that is, such an observer 
would see no very high positive free-fall frequency excitations. The reason is that the 
collapse process occurs on the much longer time scale of the Schwarzschild radius r^, so 
the modes with frequency much higher than l/r^ should remain in their ground state. 
We therefore assume that the wavepackets T"*" and T~ are in their ground states. 



where ()6.10p is used in the first and second lines, ()4.13p is used in the third line, and 
()4.1()j) is used in the last step. The problem has thus been reduced to the computation 
of the Klein-Gordon norm of the negative frequency part of the transmitted wavepacket 
T. This requires that we be more explicit about the form of the wavepacket. 

6.2.2 Norm of the negative frequency part &c thermal flux 

For definiteness we consider a spherically symmetric vacuum black hole in 3+1 dimen- 
sions, that is a Schwarzschild black hole. The Schwarzschild line element can variously 



a(T+)|^) = 



a(T-)|M/)=0. 



(6.10) 




(6.11) 
(6.12) 
(6.13) 
(6.14) 



= -(T-,T-)s„ 
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be expressed as 

ds^ = (l-'^)df-{l-'-^)-'dr^-r^{d9^ + smHdcp^) (6.15) 
r r 

= {l-y){dt^-drl)-r^{de^ + sm'^edip^) (6.16) 

= (l-—)dudv-r\de^ + sm'^ed(p'^). (6.17) 
r 

The first form is in "Schwarzschild coordinates" and = 2GM is the Schwarzschild 
radius. The second form uses the "tortoise coordinate" r^, defined by dr^, = dr /{l—rg/r) 
or = r + ln(r/rs — 1), which goes to — oo at the horizon. The third form uses the 
retarded and advanced time coordinates u = t — and v = t + r^, which are also called 
outgoing and ingoing null coordinates respectively. 

A scalar field satisfiying the Klein-Gordon equation (□ + m^)(p = can be decom- 
posed into spherical harmonics 

<p{t, r, eA) = Yl ^^^^^YimiO, 0), (6.18) 

Im 

where (f imit, r) satisfies the 1-1-1 dimensional equation 

{d^ - dl + Vim)ipim = (6.19) 

with the effective potential 

As r — > cx) the potential goes to im? . As r — >■ r^, the factor (r — r^) approaches zero 
exponentially as exp(r*/rs) with respect to r^,. Near the horizon ipim{t.,r^) therefore 
satisfies the massless wave equation, hence has the general form f{u) + g[v). 

Since the wavepacket P = '^ijnPim{t,r)Yim{0,(f)) is purely outgoing with support 
only at large radii at late times, near the horizon Pi^{t,r) must be only a function of 
the 'retarded time u — t — r*. That is, there can be no ingoing component. Since 
the metric is static, i.e. invariant with respect to f-translations, we can decompose any 
solution into components with a fixed t-frequency uj. A positive frequency outgoing 
mode at infinity has ^-dependence exp(— ia;i), hence its form near the horizon must be 
exp(— ic<;-u). 

Consider now a late time outgoing positive frequency wavepacket P that is narrowly 
peaked in Killing frequency u. Propagating backwards in time, T is the part of the 
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wavepacket that is squeezed up against the horizon, and its Yim component has the form 
Tim ~ exp{—iu!u) for all l,m. The coordinate u diverges as the horizon is approached. 
It is related to the proper time r of a free-fall observer crossing the horizon at r = 
via r ~ — To exp(— km), where n = l/2rs is the surface gravity of the black hole and the 
constant tq depends on the velocity of the free-fall observer.^ Hence the r-dependence 
of the wavepacket along the free-fall worldline is 

T ~ exp (^i^ln(-r)) (6.21) 

for r < 0, and it vanishes for r > 0. 

To find the positive frequency part we use a method introduced by Unruhj^, which 
exploits the fact that a function analytic and bounded as |r| ^ oo in the lower half 
complex r plane has purely positive frequency (see the discussion after Eqn. ()6.6p ). The 
positive frequency extension of T(r) from r < to r > is thus obtained by analytic 
continuation of ln(— r) in the lower half complex r-plane. This continuation is given by 
Inr -|- in, provided the branch cut of Inr is taken in the upper half-plane. The positive 
frequency extension of T(r) to r > is therefore obtained by replacing ln(— r) with 
Inr + 27? in ()6.21|) . which yields T(— r) exp(— ttcj/k) for r > 0. Similarly, the negative 
frequency extension of ln(— r) is given by Inr — ztt, provided the branch cut of Inr is 
taken instead in the lower half-plane. The negative frequency extension of T(r) to r > 
is therefore T(— r) exp{+TTuj / k) . Knowing these two extensions, we proceed as follows. 

Define a new wavepacket T, with support only inside the horizon, by "flipping" the 
wavepacket T{u) across the horizon (see Fig. |2)). That is, T vanishes outside the horizon 
and inside is constant on the outgoing null lines, with T(r) = T(— r) for r > 0. The 
above argument shows that the wavepackets 

T+ = c+(T + e-""/'^f) (6.22) 
T" = c_(T + e+'''^/''f) (6.23) 

have positive and negative free-fall frequency respectively. The two constants c± can 
be chosen so that T"*" + T~ agrees with T outside the horizon and vanishes (as does T) 
inside the horizon. This yields c_ = (1 — e^'^'^/'')"^ and c + /c_ = — e^'^'^/''. 

Now (T, T) = (since the two wavepackets do not overlap) and (T, T) = — (T, T) 
(since the flipped wavepacket has the reverse r-dependence), so using ()6.23|1 one finds 

(T-,T-) = ^^^^. (6.24) 



^This can be obtained by noting from the third form of the hne element (|6.17|l that along a timelike 
line (1 — rs/r)uv — 1, where the dots represent the proper time derivative. As the horizon is crossed v 
is finite, hence w ~ (r — r^)"-^ ^ e~''*/''= = e("-'')/2'-s ^ gK«_ 
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g + TCCO/K 'Y 



Figure 2: Spacetime sketch of phase contours of the transmitted wavepacket T and its flipped 
version T on either side of the horizon (dashed hue). The upper and lower signs in the exponent 
of the factor exp(=p7ra;/K) yield the positive and negative free-fall frequency extensions of T. 



Inserting this in the expression ()6.14|) for the number operator expectation value yields 
Hawking's result, 

(^1^(^)1^) = ^JJ^TT- (6-25) 

The number expectation value ()6.25|) corresponds to the result for a thermal state at 
the Hakwing temperature Th = hnj^Tx, multiplied by the so-called "greybody factor" 

r=(T,T). (6.26) 

This factor is the probability for an excitation described by the the wavepacket P to 
pile up just outside the event horizon when propagated backwards in time, rather than 
being scattered back out to infinity. Equivalently, T is the probability for the excitation 
with wavepacket P to fall across the horizon when sent in forwards in time. This 
means that the black hole would be in detailed balance with a thermal bath at the 
Hawking temperature. Yet another interpretation of T is the probability for an excitation 
originating close to the horizon with normalized wavepacket T/(T, T)^/^ to escape to 
infinity rather than scattering back and falling into the black hole. 
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6.2.3 The quantum state 

The local free-fall vacuum condition can be used to find the quantum state of 

the near horizon modes. We can do that here just "pair by pair", since the field is 
noninteracting so the state is the tensor product of the states for each outside/inside 
pair. 

Using ()6.23|) the vacuum conditions become 

a(T+)|^) oc [a(T) -e-""/V(T*)j 1^) = (6.27) 
a(T^)l^) oc [-a"f(T) + e+"'^/''a(T*)j 1^) = (6.28) 

(here we use * instead of "bar" for complex conjugation of T for typographical reasons). 
Note that since T and T have negative norm (as explained before ()6.24|l ). we have 
replaced the corresponding annihilation operators by minus the creation operators of 
their conjugates. These equations define what is called the Unruh vacuum \U) for these 
wavepacket modes. 

Now let \B) denote the quantum state of the T and T* modes such that 

a{T)\B) = (6.29) 
a{f*)\B) = 0. (6.30) 

This state is called the Boulware vacuum for these modes. In analogy with ()5.16p and 
()5.17|) the vacuum conditions ()6.27I6.28|) imply that the Unruh and Boulware vacua are 
related by 

\U) oc exp [e-"'^/''a^(f )a^(f *)j \B) (6.31) 

where the hats denote the corresponding normalized wavepackets. 

The Unruh vacuum is thus a two-mode squeezed state, analogous to the one ()5.42|) 
found at each wavevector k in the case of cosmological pair creation. There each pair 
had zero total momentum since the background was space translation invariant. In 
the present case, each pair has zero total Killing energy, since the background is time 
translation invariant. The mode T has the same positive Killing frequency as T has 
(because the Killing fiow is symmetric under the fiipping across the horizon operation 
that defines T), however its conjugate has negative Killing frequency, and therefore 
negative Killing energy. 

The Unruh vacuum is a pure state, but because of its entangled structure it becomes 
mixed when restricted to the exterior. To find that mixed state we expand the expo- 
nential in a series. Denoting by \nR^i) the level-n excitations of the modes T and T* 
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respectively, we have 

|f/)oc J]e-"'^'^/'^|ni)|n^). (6.32) 

n 

The reduced density matrix is thus 

TrL\U){U\ oc ^e-2"--/'^|nR)(nfi|, (6.33) 

n 

a thermal canonical ensemble at the Hawking temperature. 

The essence of the Hawking effect is the correlated structure of the local vacuum 
state at short distances near the horizon and its thermal character outside the horizon. 
Given this, the outgoing flux at inflnity is just a consequence of the propagation of a 
fraction F ()6.26p of each outgoing wavepacket from the horizon to infinity. 

6.3 Remarks 

In this subsection we make a large number of brief remarks about related topics that we 
have no time or space to go into deeply. Where no references are given see the sources 
listed at the end of the Introduction. 

6.3.1 Local temperature 

The Hawking temperature refers to the Killing energy, or, since the Schwarzschild Killing 
vector is normalized at infinity, to the energy defined by a static observer at infinity. A 
static observer at finite radius will perceive the thermal state to have the blueshifted tem- 
perature Tioc = Th/\^\i where |^| is the local norm of the Schwarzschild time translation 
Killing vector. At infinity this is just the Hawking temperature, whereas it diverges as 
the horizon is approached. This divergence is due to the infinite acceleration of the static 
observer at the horizon and it occurs even for an accelerated observer in the Minkowski 
vacuum of fiat spacetime (see section IB .3 .41 below) . A freely falling observer sees nothing 
divergent. 

6.3.2 Equilibrium state: Hartle-Hawking vacuum 

A black hole will be in equilibrium with an incoming thermal fiux at the Hawking 
temperature. The state that includes this incoming fiux is called the Hartle-Hawking 
vacuum. 
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6.3.3 Stimulated emission 



Suppose that the field is not in the free-fall vacuum at the horizon (jfi.lOj) . but rather 
that there are n excitations in the mode T+, so that a^{T+)a{T+)\^) = n(T+,T+)|^). 
Then instead of ()6.14|1 the expectation value of the number operator will be 

{^\N{P)\^) = n{T,T) + {n + 1){T^,T^). (6.34) 

That is, if n quanta are present to begin with in the T"*" mode, the observer at infinity 
will observe in the P mode n + 1 times the usual number of Hawking quanta, in addition 
to n times the greybody factor ()(i.2(jj) . To produce a state in which the mode is 
occupied in standard physics one would have to send in particles of enormous energy 
just before the black hole formed|S21- As explained in section I7.4.2| however, trans- 
Planckian considerations could in principle allow stimulated emission at times long after 
the black hole formed. (This has nothing to do with the standard late time stimulated 
emission of the super-radiant modes of a rotating black holejHH^ IHlj.) 

6.3.4 Unruh effect 

The argument given above for the structure of the vacuum near a black hole horizon 
applies equally well to the Minkowski vacuum near an acceleration horizon in flat space- 
time, where it is known as the Unruh effect. From a logical point of view it might be 
better to introduce the Unruh effect first, and then export it to the neighborhood of a 
black hole horizon to infer the Hawking effect. However, I chose here to go in the other 
direction 

In the Unruh effect the boost Killing field = xdt + tdx (which generates hyperbolic 
rotations) plays the role of the Schwarzschild time translation, and the corresponding 
"temperature" is h/2n. The Minkowski vacuum is the analog of the Hartle-Hawking 
equilibrium state, rather than the Unruh evaporating state. A uniformly accelerated 
observer following a hyperbolic orbit of the Killing field will perceive the Minkowski vac- 
uum as a thermal state with temperature h/2Ti\$^B\- The norm is just (x^ — t'^Y^'^, 
which is also the inverse of the acceleration of the orbit, hence the local temperature is 
the Unruh temperature Tjj = ha/2TT. As (x^ — t^)^/^ oo this temperature is redshifted 
to zero, so a Killing observer at infinity sees only the zero temperature vacuum. As the 
acceleration horizon x = ±t is approached a Killing observer sees a diverging temper- 
ature. The same temperature divergence is seen by a static observer approaching the 
horizon of a black hole in the Unruh or Hartle-Hawking states (c/. sect ion IB . 3 . 1 1 ab ove) . 
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6.3.5 Rotating black hole 

A small portion of the event horizon of a rotating black hole is indistinguishable from 
that of a Schwarzschild black hole, so the Hawking effect carries over to that case as 
well. The frequency u in ()6.25|) should be replaced by the frequency with respect to 
the horizon generating Killing field dt + ^Hd,j>, where dt and are the asymptotic time 
translation and rotation Killing vectors, and Qh is the angular velocity of the horizon. 
Thus uj is replaced by — rnQn, where m denotes the angular momentum. In effect 
there is a chemical potential mflH- See e.g. [33] for a discussion of the quantization of 
the "super-radiant" modes with u — rnQn < 0. 

6.3.6 de Sitter space 

The reasoning in the black hole case applies mutatis mutandis to de Sitter spacetime, 
where an observer is surrounded by a horizon that is locally indistinguishable from a 
black hole horizon. This leads to the temperature of deSitter spacetime [23 EE] • 

6.3.7 Higher spin fields 

The Hawking effect occurs also for higher spin fields, the only difference being (1) the 
greybody factors are different, and (2) for half-integer spin fields the Fermi distribution 
rather than the Bose distribution arises for the Hawking emission. 

6.3.8 Interacting fields 

Our discussion here exploited the free field equation of motion, but the Hawking effect 
occurs for interacting fields as well. The essence, as in the free field case, is the Un- 
ruh effect, the interacting version of which can easily be established with the help of a 
Euclidean functional integral representation of the Minkowski vacuum[36j. (This result 
was found a decade earlier, at about the same time as the original Unruh effect, via 
a theoremjnZI in the context of axiomatic quantum field theory, although the interpre- 
tation in terms of the thermal observations of uniformly accelerated observers was not 
noted until later [HH]-) The direct analog for the Hawking effect involves a Euclidean 
functional integral expression for the interacting Hartle-Hawking equilibrium state (for 
an introduction see [2111101 and references therein). 

More directly, in an asymptotically free theory one can presumably use the free field 
analysis to discover the structure of the vacuum near the horizon as in the free field 
case. The propagation of the field from that point on will involve the interactions. If 



42 



the Hawking temperature is much higher than the scale A of asymptotic freedom then 
free particles will stream away from the black hole and subsequently be "dressed" by the 
interactions and fragment into asymptotic states^T] If the Hawking temperature 
is much lower than A then it is not so clear (to me at least) how to determine what is 
emitted. 

6.3.9 Stress-energy tensor 

The Unruh and Hartle-Hawking states are "regular" on the horizon, i.e. they look like 
the Minkowski vacuum at short distances. (Recall that it is precisely the local vacuum 
property ()6.10|) or ()6.27l6.28p that determines the thermal state of the outgoing modes 
at infinity.) Hence the mean value of the stress energy tensor is finite. The Boulware 
state \B) referred to above is obtained by removing the T, T* pair excitations. This 
produces a state with a negative mean energy density that diverges as the horizon is 
approached. In fact, if even just one Hawking quantum is removed from the Unruh state 
\U) to obtain the Boulware state for that mode, a negative energy density divergence 
will be produced at the horizon. This can be viewed as the result of the infinite blueshift 
of the negative energy "hole" . That is quite odd in the context of flat space, since the 
Minkowski vacuum should be the lowest energy state, hence any other state should have 
higher energy. The explanation is that there is a positive energy density divergence on 
the horizon that more than compensates the negative energy off the horizon 

6.3.10 Back-reaction 

As previously noted the Unruh state corresponds to an entangled state of positive and 
negative Killing energy excitations. As the positive energy excitations escape to infinity, 
there must be a corresponding negative energy flux into the black hole. Studies of the 
mean value of the stress tensor confirm this. Turning on the gravitational dynamics, this 
would lead to a mass loss for the black hole via the Einstein equation. The backreaction 
driven by the mean value is called the "semi-classical" evolution. There are quantum 
fluctuations about this mean evolution on a time and length scale of the Schwarzschild 
radius, unless a large number of matter fields is invoked to justify a large limit that 
suppresses the quantum fluctuations. 

6.3.11 Statistical entropy 

The entangled structure ()fi.32j) of the Unruh state leads to a mixed state ()(i.33|l when ob- 
servations are restricted to the region outside the horizon. The "entanglement entropy" 
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—Trplnp of this state is the same as the thermal entropy of the canonical ensemble 
()6.33j) . Summing over all modes this entropy diverges due to the infinite density of 
modes. To characterize the divergence one can use the thermodynamic entropy density 
s oc of a bath of radiation at temperature T. The local temperature measured by a 
static observer (c/. section I?). 3. 1|1 is given by 



where i is the proper distance to the horizon on a surface of constant t, and the relation 
K = {d\^\/di)H has been used. Thus the entropy diverges like 



where ic is a cutoff length above the horizon.® 

What is the meaning of this entropy? It seems clear on the one hand that it must 
be included in the black hole entropy, but on the other hand it must somehow be mean- 
ingfully cut off. It is natural to try to understand the scaling of the black hole entropy 
with area in this way, but this is only the contribution from one quantum field, and 
there is also the classical contribution from the gravitational field itself that emerges 
from the partition function for quantum gravity [Hj. The apparent dependence on the 
number of fields is perhaps removed by the corresponding renormalization of Newton's 
constant jlH]. (The last reference in |3H] is a review.) However this is regularization 
dependent and hence difficult to interpret physically, and moreover at least in dimen- 
sional regularization vector fields and some non-minimally coupled scalars contribute 
negatively to renormalizing G, which does not seem to match the entanglement entropy. 
There has been much work in this area but it remains to be fully understood. 

6.3.12 Information loss 

Two types of potential information loss occur in black hole physics. First, when some- 
thing falls into a black hole any information it carries is apparently lost to the outside 
world. Second, when a black hole radiates Hawking quanta, each radiated quantum is 
entangled with a partner lying inside the horizon, as Eqn. ()6.32|) shows. As long as the 

^Sorkin*^ introduced the notion of black hole entanglement entropy, and with collaborators 
computed it in the presence of a regulator. A mode-by-mode version of the thermal entropy calculation 
was first done by 't Hooft^Hlj called it the "brick wall model" because of a Dirichlet boundary 
condition applied at ic- 



Tioc = Th/\^\ ^ Th/kE = l/2ni, 



(6.35) 




(6.36) 



44 



black hole does not evaporate completely, all information remains available on a space- 
like surface that crosses the horizon and enters the black hole. If on the other hand the 
black hole evaporates completely, then no single spacelike surface stretching to infinity 
and filling "all space" can capture all information, according to the semi-classical anal- 
ysis of the Hawking effect. Rather a disconnected surface behind the horizon must be 
included. The information on this disconnected surface flows into the strong curvature 
region at the singularity, where its fate is not yet understood. 

This situation has generated much discussion. Some researchers (myself included) 
see the information loss to the outside world not as a sign of the breakdown of quantum 
mechanics but just as a consequence of the mutabihty of spatial topology in quantum 
gravity. When a black hole is about to evaporate completely it looks very small to the 
world outside the horizon. However this outside smallness has absolutely nothing to do 
with the size of the region inside available for storing information. Let us look into this 
a bit more. 

Consider for example the spacelike singularity at r = inside the Schwarzschild black 
hole. The metric in Eddington-Finkelstein coordinates is 

ds'^ = (1 - rs/r)dv'^ - 2dvdr - r'^{de^ + sin^ 9d^'^), (6.37) 

where v = t + is the advanced time coordinate defined below ()6.17p . Inside the 
horizon, where r < r^, a line of constant r, 6, is spacelike, and has a proper length 
L(r, Av) = {vs/r — lY^'^Av. This goes to infinity as the singularity is approached, for 
any interval of advanced time Av. Hence there is no dearth of space inside. On the 
other hand, the transverse angular dimensions go to zero size, and also we do not know 
how to describe the spacetime too near to the singularity. Therefore let's stop at the 
radius where the curvature ~ r^/r^ is equal to the Planck curvature, i.e at r ~ rl^^ in 
Planck units. Then then proper length goes as L{ry^,Av) ~ rl^^Av. For a solar mass 
black hole we have ~ 1km ~ 10^^ in Planck units, so rl^^ ~ 10^^. That means that 
for an external advanced time of Av = 1 second, the proper length inside is one million 
light years. After a day or so, the length is the size of the universe, and so on over the 
Hawking lifetime ~ 10^^^. 

1/3 

What happens to the future of this spacelike r = rs cylinder inside the black hole 
is governed by quantum gravity. We don't know what form the evolution takes. It is 
conceivable that time just stops running, like a frozen engine, producing a boundary of 
spacetime. It seems much more likely however that spacetime persists beyond, either 
into "quantum foam" or into a plump baby universe (see e.g. ^H] and references therein) 
or universes |5Uj. In any of these scenarios, the tiny outside size of a black hole at the 
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end stage of Hawking evaporation is no indication of the information carrying capacity 
of the interior. 

Others beheve that the vahdity of quantum mechanics requires that all information 
winds up available on the exterior slice after the evaporation. Results from string theory 
are often invoked to support this viewpoint. (For a critique of these arguments see jKT].) 
If true, it would require some breakdown of the semi-classical description where none 
seems to be otherwise called for. The role of ultra high frequencies in the Hawking effect 
is often brought up in this context, however they are irrelevant since the derivation of 
the correlated structure of the vacuum and the Hawking effect does not need to access 
those frequencies. For more on this see the discussion of the trans-Planckian question 
in section [7| 

6.3.13 Role of the black hole collapse 

The Hawking flux continues in a steady state long after the collapse that forms a black 
hole, which suggests that the collapse phase has nothing to do with the Hawking effect. 
Indeed the derivation given above makes use of only the free-fall vacuum conditions near 
the horizon, and the collapse phase plays no role. On the other hand if, as Hawking orig- 
inally did, we follow the positive Killing frequency wavepacket T all the way backwards 
in time, it would go through the collapse phase and back out to infinity where Killing 
frequency and free-fall frequency coincide. Were it not for the time dependence of the 
background during the collapse, there would be no change of the Killing frequency, so 
there would be no negative frequency part of the ingoing wavepacket and there would 
be no particle creation. As discussed in the following section this propagation through 
the collapse phase invokes ultra high frequency field modes, and therefore may not even 
be physically relevant. However it seems clear is that whatever physics delivers the out- 
going vacuum near the horizon, it must involve some violation of the time translation 
symmetry of the classical black hole background, even if it does not involve the collapse 
phase. In the lattice model of section ESI the violation arises from microscopic time de- 
pendence of the lattice spacing. In quantum gravity it may come just from the quantum 
gravitational fluctuations (see section 17.4. 3|) . 

7 The trans-Planckian question 

A deep question arises on account of the infinite redshift at the black hole horizon: do 
the outgoing modes that carry the Hawking radiation really emerge from a reservoir 
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of modes with frequency arbitrarily far beyond the Planck frequency just outside the 
horizon, or is there another possibility? Does the existence or properties of the Hawking 
effect depend on the existence of such a trans- Planckian reservoir! 

The reasoning leading to the expression ()6.32j) for the Unruh state in terms of positive 
Killing frequency modes is largely shielded from this question. The essential input is the 
free-fall vacuum conditions ()6.10|) . which can be applied at any length or time scale much 
shorter than the Schwarzschild radius and inverse surface gravity (which are roughly the 
same unless the black hole is near extremally rotating or charged). There is no need to 
appeal to Planckian or trans-Planckian frequencies. 

From this persepctive it is clear that, as far as the derivation of the Hawking effect 
is concerned, the only question is whether or not the free-fall vacuum in fact arises at 
short distances near the horizon from the initial conditions before collapse. As mentioned 
earlier, the modes with frequencies much higher than the inverse of the collapse time 
scale would be expected to remain unexcited. Nevertheless, in standard relativistic field 
theory these modes arise from trans-Planckian modes. 

Consider an outgoing wavepacket near the horizon and peaked around frequency uf 
as measured by a free-fall observer crossing the horizon at the advanced time vi. At an 
earlier time f 2 = f i — Av, the wavepacket would be blueshifted and squeezed closer to 
the horizon, with exponentially higher free-fall frequency. 

For a solar mass black hole and Av = 2 seconds, the ratio is exp(lO^). 

To predict the state of the positive free-fall frequency modes and from the 
initial state thus seems to require trans-Planckian physics. This is a breakdown of the 
usual separation of scales invoked in the application of effective field theory and it leaves 
some room for doubt jHS EOl EB] about the existence of the Hawking effect. 

While the physical arguments for the Hawking effect do seem quite plausible, the 
trans-Planckian question is nevertheless pressing. After all, there are reasons to suspect 
that the trans-Planckian modes do not even exist. They imply an infinite contribu- 
tion to black hole entanglement entropy from quantum fields, and they produce other 
divergences in quantum field theory that are not desirable in a fundamental theory. 

The trans-Planckian question is really two-fold: 

1. Is the Hawking effect universal, i.e. insensitive to short distance physics, or at 
least can it be reliably derived in a quantum gravity theory with acceptable short 
distance behavior? 
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2. If there is no trans-Planckian reservoir, from where do the outgoing black hole 
modes arise? 

7.1 String theory viewpoint 

String theory has made impressive progress towards answering the first question, at 
least for some special black holes. In particular |54t 133] . some near-extremal black holes 
become well-understood D-brane configurations in the weak coupling limit, and super- 
symmetry links the weak coupling to strong coupling results. Thus the Hawking effect 
can be reliably analyzed in a full quantum gravity theory. By a rather remarkable and 
unexpected correspondence the computations yield agreement with the semi-classical 
predictions, at least in the long wavelength limit. Moreover, the D-brane entropy is 
understood in terms of the counting of microstates, and agrees with the correspond- 
ing black hole entropy at strong coupling, just as the supersymmetry reasoning says 
it should. From yet another angle, the AdS/CFT duality in string theory offers other 
support IHEl • There the Hawking effect and black hole entropy are interpreted in terms of 
a thermal state of the CFT (a conformally invariant super- Yang- Mills theory). However, 
neither of these approaches from string theory has so far been exploited to address the 
origin of the outgoing modes, since a local spacetime picture of the black hole horizon 
is lacking. This seems to be a question worth pursuing. 

7.2 Condensed matter analogy 

Condensed matter physics provides an analogy for effective field theory with a fundamen- 
tal cutoff, hence it can be used to explore the consequences of a missing trans-Planckian 
reservoir. (For a review of these ideas see [37], and for a very brief summary see [3H].) 
The first such black hole analog was Unruh's sonic black hole, which consists of a fluid 
with an inhomogeneous flow exceeding the speed of sound at a sonic horizon. A molec- 
ular fluid does not support wavelengths shorter than the intermolecular spacing, hence 
the sonic horizon has no "trans-molecular" reservoir of outgoing modes. Unruh found 
that nevertheless outgoing modes are produced, by a process of "mode conversion" from 
ingoing to outgoing modes. This phenomenon comes about already because of the alter- 
ation of the dispersion relation for the sound waves. It has been studied in various field 
theoretic models, however none are fully satisfactory since the unphysical short distance 
behavior of the field is always eventually called into play. The model that most closely 
mirrors the fluid analogy is a falling lattice jSH] which has sensible short distance physics. 
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The mode conversion on the lattice involves what is known as a Bloch oscillation in the 
condensed matter context. Here I will briefly explain how it works. 



7.3 Hawking effect on a falling lattice 



The model is 2d field theory on a lattice of points falling freely into a black hole. We 
begin with the line element in Gaussian normal coordinates, 



dt^ — a^{t, z) dz^ . 



(7.2) 



A line of constant z is an infalling geodesic, at rest at infinity, and the "local scale factor'' 
a(t, z) satisfies 



a(t, z oo) 
a(0, z) 
a{t,ZH{t)) 



1 
1 

K.t 



for Kt > 1, 



(7.3) 
(7.4) 
(7.5) 



expressing the facts that the metric is asymptotically fiat, the coordinate z measures 
proper distance on the t = time slice, and at the horizon znit) the time scale for 
variations of the local scale factor is the surface gravity k. The specific form of a{t, z) 
is not required for the present discussion. (For details see ^59J.) A scalar field on this 
spacetime is governed by the action 



S 



^ / d'^x 

1 
2 



dtdz 



a{z,t){dtipf 



a{t, z] 



Now we discretize the z coordinate with spacing 6: 



m6. 



d^(p Dip 



5 



(7.6) 
(7.7) 

(7.8) 



where ipm = f{zm)- The proper distance between the lattice points Zm and Zm+i on a 
constant t slice is approximately a{t, Zm)S. At t = this is just 6 everywhere, that is the 
points start out equidistant. However, since they are on free fall trajectories at different 
distances from the horizon, they spread out as time goes on. In particular the lattice 
spacing at the horizon grows with time like ~ Kt. 
The discretized action is 



^1 



lattice 



dtj2 



am{t){dtiPr. 



(t) + a^{t)\ 



(7.9) 
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The discrete field equations produce the dispersion relation 

u^ik) = ±— ^ sin(M/2) (7.10) 
a{z, t)o 

for a mode of the form exp(— icu^t + iA;z), provided dta <C and dza -C k. (See Fig. IHl) 
For small wave numbers the lattice dispersion agrees with the continuum, however it is 




Figure 3: Dispersion relation ujS = ±2 sin(/c(5/2) plotted vs. k5. Wavevectors differing by 
271/6 are equivalent. Only the Brillouin zone \k\ < 11/6 is shown. 



periodic in translation of k by 211/6. At the wavenumber k = n/S there is a maximum 
frequency 2/a6 and vanishing group velocity duo/dk . Beyond that wavenumber the 
group velocity reverses, and k is equivalent to k — 271/5 which lies in the Brillouin zone 
\k\ < 7x/b. 

On the lattice the trans-Planckian redshift cannot take place, because of the lat- 
tice cutoff. Hence the outgoing modes — provided they exist — must come from ingoing 
modes. Both WKB "eikonal" trajectory and numerical evolution of the discrete wave 
equation confirm that indeed this occurs. The behavior of a typical wavepacket through- 
out the process of bouncing off the horizon is illustrated in Fig. 0] The real part of the 
wavepacket is plotted vs. the static coordinate at several different times. Following 
backwards in time, the packet starts to squeeze up against the horizon and then a trail- 
ing dip freezes and develops oscillations that grow until they balloon out, forming into 
a compact high frequency wavepacket that propagates neatly away from the horizon 
backwards in time. 

The mode conversion can be understood as follows, following a wavepacket peaked 
around a long wavelength A ^ 5 backwards in time. The wavepacket blueshifts as it 
approaches the horizon, eventually enough for the lattice structure and therefore the 
curvature of the dispersion relation to be felt. At this point its group velocity begins 
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Figure 4: A typical wavepacket evolution on the lattice. The oscillations of the incoming 
wavepacket are too dense to resolve in the plots. 

to drop. In the WKB calculation the wavepacket motion reverses direction at a turning 
point outside the horizon. This occurs before its group velocity in the falling lattice 
frame is negative, so it is faUing in at that stage because its outward velocity is not 
great enough to overcome the infaUing of the lattice. 

As the wavepapcket continues backwards in time now away from the black hole its 
wavevector continues to grow until it goes past the edge of the Brillouin zone, thus 
becoming an ingoing mode also in the lattice frame. This reversal of group velocity is 
precisely what happens in a "Bloch oscillation" when a quantum particle in a periodic 
potential is accelerated. As the turnaround is occuring, another equally important ef- 
fect is that the time-dependence of the underlying lattice is felt. Thus, unlike in the 
continuum limit of this stationary background, the Killing frequency of the wavepacket 
is no longer conserved. 

Following the wavepacket all the way backwards in time out to the asymptotic region, 
it winds up with a short wavelength of order S and a large frequency of order 1/S. This 
frequency shift is absolutely critical to the existence of the outgoing modes, since an 
ingoing low frequency mode would simply fall across the horizon. Only the "exotic" 
modes with sufficiently high frequency will undergo the mode conversion process. 

The Hawking flux is determined by the negative frequency part of the ingoing 
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wavepacket. The eikonal approximation just described does not capture the negative 
frequency mixing which occurs during the turnaround at the horizon. Just as in the 
continuum, the wavepacket squeezed against the horizon has both positive and negative 
free-fall frequency components. As these components propagate backwards in time away 
from the horizon, their frequency slowly shifts, but their relative amplitude remains fixed. 
Hence the norm of the negative freqency part of the incoming wavepacket turns out to 
be just what the continuum Hawking effect indicates, with small lattice corrections. Put 
differently, the infalling vacuum is adiabatically modified by the underlying microscopic 
time dependence of the lattice in such a way that the Unruh conditions ()6.10|) on the 
state of the outgoing modes at the horizon are satisfied. 

7.4 Remarks 

7.4.1 Finite Entanglement entropy 

Since the lattice has a short distance cutoff the entanglement entropy between the modes 
just inside and outside of the horizon is finite at any time. As illustrated in Figure any 
given entangled pair of vacuum modes began in the past outside the horizon, propagated 
towards the horizon where it was "split", and then separated, with one half falling in 
and the other converted into an outgoing mode. As time goes on, new pairs continually 
propagate in and maintain a constant entanglement entropy. 



Figure 5: The ancestors of a Hawking quantum and its negative energy partner. In stan- 
dard relativistic field theory the ancestors are trans-Planckian and pass through the collapsing 
matter at the moment of horizon formation. On the lattice the ancestors ar Planckian and 
propagate in towards the black hole at late times. 




planckian 
\ ancestors 
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7.4.2 Stimulated emission of Hawking radiation at late times 

As discussed in section 16 .3 .3^ if the Unruh vacuum conditions fl(j.l()|l do not liold at 
the horizon then stimulated emission of Hawking radiation will occur. In the falling 
lattice model, these horizon modes arise from ingoing modes long after the black hole 
formed as shown in Figure Thus it is possible to stimulate the emission of Hawking 
radiation by sending in radiation at late times, in contrast to the usual continuum 
case of a static black hole. This seems a generic feature of theories with a cutoff, for 
which the outgoing modes must arise from modes that are ingoing after the collapse. 
(Something like it should happen also in string theory if, as many suppose, the trans- 
Planckian reservoir at the horizon is also eliminated there.) Note however that the linear 
model described here is surely a gross oversimplification. Turning on the gravitational 
interactions between the modes and the background, one is led to a picture in which the 
modes "dissipate" when propagated backwards in time to the Planckian regime. Hence 
what really produces the outgoing Hawking quantum must be a complicated collective 
mode of the interacting vacuum that "anti-dissipates" as it approaches the horizon and 
turns around. Calculations exploring this in quantum gravity were carried out in [60] . 

7.4.3 Lattice time dependence and geometry fluctuations 

The microscopic time dependence of the lattice, i.e. the slow spreading of the lattice 
points, plays a critical role in transforming an ingoing mode with high Killing frequency 
to an outgoing mode with low Killing frequency, and in allowing for mixing of positive 
and negative frequencies despite the stationarity of the continuum black hole back- 
ground. This suggests the conjecture that in quantum gravity the underlying quantum 
fiuctuations of the geometry that do not share the stationarity of the black hole back- 
ground metric might play this role. A step towards understanding this might be provided 
by a two or three dimensional version of the lattice model, in which the density of lattice 
points remains fixed but their microscopic positions fiuctuate. This is precisely what 
happens with an inhomogeneous fiow of a real molecular fiuid. The quantum gravity 
analysis of [^0] lends some support to this conjecture, though it is not clear whether 
the Lorentz violation seen there is just an artifact of a non-covariant cutoff or a feature 
introduced by the global geometry of the black hole spacetime. 
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7.5 Trans-Planckian question in cosmology 

The modes producing the inflationary perturbation spectrum (c/. section IB. 4. HI redshift 
exponentially from their trans-Planckian origins. It has been suggested that this might 
leave a visible imprint on the perturbation spectrum, via a modified high frequency 
dispersion relation and/or a modified initial quantum state for the field modes. As il- 
lustrated by the Hawking process on the lattice however, as long as the redshifting is 
adiabatic on the timescale of the modes, they would remain in their ground states. If 
the Hubble rate H during inflation were much less than the Planck mass Mp (or what- 
ever scale the modified dispersion sets in) the modes could be treated in the standard 
relativistic manner by the time the perturbation spectrum is determined. At most one 
might expect an effect of order H/Mp, and some analyses suggest the effect will be even 
smaller. This all depends on what state the modes are in. They cannot be too far from 
the vacuum, since otherwise running them backwards in time they would develop an 
exponentially growing energy density which would be incompatible with the inflationary 
dynamics. Hence it seems the best one can say at present is that there may be room for 
noticable deviations from the inflationary predictions, if H/Mp is large enough. (See for 
example [HI] and references therein for discussions of these issues.) 
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